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ABSTRACT 

Context. Understanding the collisional properties of ice is important for understanding both the early stages of planet formation and 
the evolution of planetary ring systems. Simple chemicals such as methanol and formic acid are known to be present in cold protostellar 
regions alongside the dominant water ice; they are also likely to be incorporated into planets which form in protoplanetary disks, and 
planetary ring systems. However, the effect of the chemical composition of the ice on its collisional properties has not yet been studied. 
Aims. Collisions of 1.5 cm ice spheres composed of pure crystalline water ice, water with 5% methanol, and water with 5% formic 
acid were investigated to determine the effect of the ice composition on the collisional outcomes. 

Methods. The collisions were conducted in a dedicated experimental instrument, operated under microgravity conditions, at relative 
particle impact velocities between 0.01 and 0.19 m s _1 , temperatures between 131 and 160 K and a pressure of around 10~ 5 mbar. 
Results. A range of coefficients of restitution were found, with no correlation between this and the chemical composition, relative 
impact velocity, or temperature. 

Conclusions. We conclude that the chemical composition of the ice (at the level of 95% water ice and 5% methanol or formic acid) 
does not affect the collisional properties at these temperatures and pressures due to the inability of surface wetting to take place. At a 
level of 5% methanol or formic acid, the structure is likely to be dominated by crystalline water ice, leading to no change in collisional 
properties. The surface roughness of the particles is the dominant factor in explaining the range of coefficients of restitution. 

Key words, accretion, accretion disks - astrochemistry - planets and satellites: formation - planets and satellites: rings - protoplan¬ 
etary disks 


1. Introduction 

Water ice is abundant in the interstellar medium and has also 
been observed towards dense clouds, cloud cores, protostellar 
regions and protoplanetary disks ( Oberg et al.|201 la| . Recentre- 
search studying deuterium-to-hydrogen enrichments in the solar 
system has shown that water ice was present in the solar neb¬ 
ula protoplanetary disk (Cleeves et al. 2014 1 . Water ice is also 
present in planetary ring systems such as the rings of Saturn 
|Cuzzi & Pollack] 1978| |Cuzzi et al. [ 1980] >, where its collisional 
properties play an important role in the structure and dynamical 
evolution of the rings. It is not clear whether these rings are old 
or young; in the former case, the rings will be residues from the 
formation era. The presence of water ice across star and planet 
forming regions may therefore influence the processes of planet 
formation. 

It has been proposed that planets form from the dust in a 
protoplanetary disk by a process of dust aggregation ( |Weiden-| 
schilling} 1977 19801. While this theory has gained widespread 


(the so-called bouncing barrier) ( |Guttler et al.|2010| |Zsom et aT] 
2010| Kothe et al.|2013|l. In recent years, the presence of ice has 


is a critical velocity, reducing with particle size (Weidenschilling 
|1977]>, beyond which particles tend to bounce rather than stick 


acceptance, there remains a problem in that growth between cen¬ 
timetre and kilometre sizes has not been demonstrated, as there 


been postulated as one possible solution to the bouncing bar¬ 
rier; it has been shown that ice particles have a larger adhesive 
and rolling friction force ( Gundlach et al.|201 1 1, a higher stick¬ 
ing thresh old (|Gundlach & Blum||2014| l and reduced elasticity 
( |Hertz sch 2002J compared to silicate dust particles which will 
increase the threshold velocity for sticking. In addition, recent 
model simulations have shown that ice condensation could en¬ 
able dust grains to grow to decimetre sizes around the snowline 
( Ros & Johansen|2013 i and even to icy planetesimals if the ini¬ 
tial ice grains were submicrometre-sized (Kata oka et al.|20f3) . 

In the context of this paper, ring systems are also of inter¬ 
est. Until recently, rings had exclusively been known to exist 
around the four outermost solar system planets, though with the 
detection of two small dense rings around the Centaur (10199) 
Chariklo ( |Braga-Ribas et al.|2014[|Duffard et al.|2014| >, it seems 
that ring systems are not exclusively present around giant plan¬ 
ets. The rings around Saturn are the most well studied and con¬ 
sist predominantly of small icy bodies (Cuzzi et al. [2010|, be¬ 
tween 1 cm and 10 m diameter, dominated by water ice (|Zebker 
et a!.11985| >. The ring dynamics are complex, influenced by per¬ 


turbations by resonances between the ring particles and nearby 
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moons, which are counteracted by frequent inelastic collisions 
at low relative velocities. Typical collision velocities in the un¬ 
perturbed rings are well below 1 cm s -1 (Esposito |2002 


Col- 


|well e t al. 2009), but gravitational perturbations may raise the 
average collision speed considerably. The collisions dissipate ki¬ 
netic energy from the rings and thereby determine the stability of 
the rings ( Salo et al.| 2001 | |Schmidt et al.|| 200 T| ); consequently, 
data pertaining to collisional properties and outcomes are vital 
to model efforts which are used to understand ring instabilities, 
wakes and overall structure, including the very strongly confined 
heights observed in Saturnian and other ring systems (Colwell 


|et al.|2006| [20071 |Hedman et al.|2007[ [Thomson et al.|2007). 

While water is abundant in planet forming regions, many 
other molecules have also been detected in protostellar regions. 
Consequently, the icy mantles on dust grains in protoplanetary 
disks are likely to be composed of water dominated ices which 
may also contain a plethora of "contaminants", including many 
simple molecules such as CO, CO 2 , CH 3 OH, CH 4 , NH 3 and 
HCOOH (e.g. Pontoppidan et al.] ( 2004| ), [Aikawa et al.| ( |2012] >, 


tution was thought to be due to the difference in velocity ranges 
between the two studies. 

The range of coefficients of restitution and lack of depen¬ 
dence on impact velocity and normalised impact parameter is in 
contrast to previous studies of ice collisions which show a de¬ 
crease in coefficient of restitution with increasing impact veloc¬ 
ity ( |Bridges et al.| ( T984| ), |Hatzes et al.| (fl988)) and very little loss 
of kinetic energy for glancing collisions (i.e. b/R ~ 1) (Supulver 
et al. 1995| . The reason for this was thought to be the differ¬ 
ence in surface roughness; while the ice surfaces in the work of 
Bridges et al. (1984]), Hatzes et al. ( 1988|) and Supulver et al. 


(|1995|l were smooth, the ice surfaces in the studies of HeiBel 


Obe rg et aI7| ( |201 lb] f Noble et al. ( |2013 1 ). Recent work has sug 


gested that the most abundant volatile species in planets of solar 
composition are H2O, CO, CO 2 , CH3OH and NH3 (Marboeuf| 
|et al.|2014| . ft is also possible that some of these species may be 
present in planetary rings such as those of Saturn. While Saturn’s 
rings are composed of mainly water ice, there is evidence that 
other species are present ( |Poulet et al.|2003j |Cuzzi et al.[2 009) 
and it is not unreasonable to suppose that some of those may in¬ 
clude the above molecules. This provides the main motivations 
for this paper: firstly to study the effects of minor chemical pol¬ 
lutants on collisions between icy particles at low velocities, and 
secondly to understand how the outcomes of these collisional 
processes affect processes of planet formation and ring system 
evolution. 

The work on this paper builds on our previous publication 


|mann et ah] (2010} and Paper I were rough and anisotropic. This 
explains the ranges of coefficient of restitution and the lack of 
dependence on impact velocity and normalised impact parame¬ 
ter and is in accordance with the work of |Hatzes et al.| ( fi988j l 
where both frost and roughened surfaces were found to reduce 
the coefficient of restitution. For a full comparison of the exper¬ 
imental conditions in these studies, we refer the reader to Paper 
I. 

An inherent limitation in the studies of ice collisions to date 
(as related to planet formation and planetary ring systems) is 
their focus on pure water ice. It is possible that the presence of 
other chemicals might change the collisional properties; [Bridges 


|et al.| ( [T996] > studied the sticking forces between frost coated wa¬ 
ter ice at atmospheric pressure and temperatures between 110 - 
150 K and discovered that methanol frosts had stronger sticking 
forces than water frosts. The melting point of methanol is 176 K 
at atmospheric pressure so it is possible that the increased stick¬ 
ing forces are due to surface melting. This is the same mecha¬ 
nism that leads to ice aggregation in atmospheric clouds (Hobbs] 
1965[), leading to increased aggregation around the melting point 


of water ( |Hobbs et al.|1974| l. In protoplanetary disks and plan¬ 
etary rings, the temperatures and pressures are too low for this 
process to dominate. Both methanol and formic acid have been 


of millimetre-sized ice particles (both spheres and irregularly 

detected in protostellar regions in the solid phase ( 

Pontoppidan 

shaped fragments, diameters 4.7-10.8 mm) were studied under 

et al. 

(2004), (Oberg et al.|(2011b|), Aikawa et al. 

u 20 l 2 ), Schutte 


microgravity conditions on a parabolic flight. The main param¬ 
eters studied were the coefficient of restitution, e, and the im¬ 
pact parameter, b. The coefficient of restitution is the ratio of the 
relative velocities of the particles after and before the collision, 
which is related to the translational kinetic energy lost in the col¬ 
lision: 


et al. 


Va 

e — — 
Vb 


(1) 


where a and b denote after and before the collision respectively. 
The impact parameter is the distance of closest approach of 
the two particles perpendicular to their relative velocity vector. 
Henceforth we use the normalised impact parameter b/R , which 
is the impact parameter, b, normalised to R, the distance between 
the centre of masses of the two particles at the point of collision, 
i.e. the sum of the two radii for spherical particles. 

The relative collision velocities ranged from 0.26 to 
0.51 m s" 1 . The coefficients of restitution were spread evenly 
between 0.08 and 0.65 with no dependence of this property 
on either impact velocity or normalised impact parameter. The 
spread of coefficients of restitution was attributed to the surface 
roughness of the particles. This corroborates previous work by 
HeiBelmann et al. (2010 1 , where collisions of 1.5 cm (diameter) 
ice spheres were studied under microgravity conditions. In this 
case, the impact velocities were between 0.06 and 0.22 m s _1 and 
the coefficient of restitution was spread evenly between 0.06 and 
0.84. The difference between the ranges of coefficients of resti¬ 


(1999[ ), |Zasowski et al.| (20091 )) at the 1-30% abundance 
level relative to water ice for methanol and 1-5% level for formic 
acid (Boogert e t ak]|2008| . We have therefore chosen to focus 
this study on icy particles formed by freezing either pure water 
or 5% methanol or formic acid in water, to investigate the effects 
of the chemical composition of ice on its collisional properties. 
While the eventual aim is to investigate the effects of the chem¬ 
ical composition on planet formation, the form of ice used here 
(solid ice spheres) is more relevant to planetary ring systems than 
protoplanetary disks, where the ice is likely to be in the form of 
porous agglomerates (e.g. |Kataoka et al.| ( |2013[ )). Therefore the 
results will be more immediately applicable to planetary rings, 
but it is hoped that they will be a starting point for discussion of 
ice in protoplanetary disks as well. 


2. Experimental details 

The experimental setup has previously been discussed in detail 
( [Salter et al.| ( j2009j l and Paper I) and is summarised again here 
for completeness (see Fig. [TJ. The experiment was designed and 
built to carry out low velocity particle collisions on parabolic 
flights. Prior to take off, the particles are loaded into a pre¬ 
cooled (to around 77 K) rotating double helix particle reservoir, 
or colosseum (1). This is rotated (2) during the flight to line up 
two diametrically opposed ice particles with motorised pistons 
(3) which accelerate the particles to a pre-defined constant ve¬ 
locity directed towards the centre of the collision volume, where 
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Fig. 1 . Computer aided design schematic of the experimental setup 
(adapted from |Salter et al. ^2009} . The particles are stored in an alu¬ 
minium colosseum (1) which is rotated (2) to line up two colosseum 
holes with two diametrically opposed motorised pistons (3) which col¬ 
lides the particles either with a 52 mm diameter ice target (4) placed at 
the centre of the collision volume or with each other when the target 
is removed. The system is passively cooled by passing liquid nitrogen 
through a cooling ring (5) prior to take off. The particles are kept within 
the colosseum and protected from radiative heating by an aluminium 
shield (6), and a 45 kg copper block (7) acts as a thermal reservoir, 
keeping the system cool during the flight. The entire set up is enclosed 
within a vacuum chamber (8). 


they collide in perfect free fall conditions either with a target (4) 
or with each other. The system is cryogenically cooled to tem¬ 
peratures around 77 K by passing liquid nitrogen through a cool¬ 
ing ring (5) prior to take off and is kept cool during flight with 
the use of a U-shaped aluminium heat shield (6) and a 45 kg 
copper block (7) on which the particle reservoir sits. The copper 
acts as both a thermal reservoir and a cryo-pump for any residual 
gas. The entire set up is encased inside a vacuum chamber (8), 
with a typical residual gas pressure of 10 s mbar. The pistons 
enter the collision volume through two differentially pumped 
feedthroughs. A collision is initiated by accelerating the pistons 
(which are controlled by Labview software) to a maximum ve¬ 
locity of between 0.01 and 0.1 m s _1 . The pistons make contact 
with the particles, are brought to a sharp stop a few millimetres 
later, and then are immediately retracted to their starting posi¬ 
tions. The piston tips remain within the cooled region at rest, 
meaning that they are at the same temperature as the particles 
and do not induce particle heating or surface melting during col¬ 
lision. A high speed camera with a frame rate of 107 frames per 
second was combined with prism optics to capture two views of 
the collision 48.8° apart. 

The particles used in this experiment were ice spheres of 1.5 
cm diameter. Some were composed of pure, distilled water, and 
others were composed of a 5% solution of methanol or formic 
acid in distilled water. The particles were produced using spher¬ 
ical moulds and a standard kitchen freezer. Making ice parti¬ 
cles in this way produces crystalline hexagonal ice with a rough, 
anisotropic surface ( Heifielmann et al.|2010) . The particles were 
removed from the freezer and placed directly into a container 
of liquid nitrogen situated within the chamber. Once the parti¬ 
cles had reached 77 K, they were individually loaded into the 
pre-cooled sample reservoir. Surface frosting was minimised by 
pre-cooling the particle reservoir to ~77 K prior to loading and 
evacuating the chamber to a base pressure of around 1CU 5 mbar 
immediately after loading while the ice particles were still out- 
gassing N 2 , thereby preventing frost from forming on the particle 


surfaces. The target was composed of pure crystalline water ice 
and was placed in the centre of the chamber immediately prior to 
evacuation. The target was placed in a different orientation with 
respect to the chamber for each flight, giving access to different 
normalised impact parameters. The normalised impact parame¬ 
ter is related to the angle of the target to the direction of travel 
(6) by the following equation: 

b/R = cos{6) (2) 


On the first flight of this campaign, the target was at 90° to the 
direction of travel of the particles, giving b/R — 0 - completely 
head on collisions. The second flight had the target at 60° to the 
direction of travel, giving b/R = 0.5. The third flight had the tar¬ 


get mounted at 30° to the direction of travel, giving b/R = V3/2. 
Technical problems prevented collection of data in this configu¬ 
ration and so particle-particle collisions were performed instead. 
It was therefore necessary to open the chamber during the flight 
to remove the target. The particles remained within the sam¬ 
ple reservoir and hence their exposure to the atmosphere of the 
plane was minimal. It is unlikely that significant frost formation 
could occur on the particle surfaces under these conditions and 
comparison with previous particle-particle collisions conducted 
under similar conditions (but without the opening of the cham¬ 
ber) shows that this assumption is correct (see Section |4~T| for a 
comparison of our data with that of HeiBelmann et al,|^2010| l). 
Fig.[2]shows an example collision of an ice sphere containing 5% 
formic acid with the ice target at 90° to the direction of travel. 

Table|T]shows a break down of all the collisions successfully 
performed during this campaign. The numbers of collisions suit¬ 
able for analysis and unsuitable for analysis are given. Residual 
acceleration of the plane (when the quality of microgravity is 
poor and the camera moves horizontally with respect to the par¬ 
ticles as a result) was a particular problem on the second day of 
the campaign (target at 60°) and hence only 8 collisions suitable 
for analysis were obtained. Poor quality hampered image analy¬ 
sis efforts in a few cases and there was one case of fragmentation. 
This is thought to be related to the way in which the ice spheres 
are formed rather than a real collisional outcome. Table [2] shows 
a break down of the collisions suitable for analysis, the results of 


which are presented in Section 4.1 


Table 1 . A break down of the collisional outcomes observed in this 
study. The angles refer to the angle between the target and the direction 
of travel of the particles. The collisions that were suitable for analysis 
are detailed in Section |4~T1 The other collisions were excluded from 
the analysis due to residual acceleration (poor quality of microgravity), 
poor image quality and fragmentation. 



Collisional outcome 

Number of collisions 

Target at 90° 

Suitable for analysis 

35 


Residual acceleration 

1 


Poor quality 

4 


Fragmentation 

1 

Target at 60° 

Suitable for analysis 

8 


Residual acceleration 

41 


Poor quality 

5 


Fragmentation 

0 

No target 

Suitable for analysis 

15 


Residual acceleration 

3 


Poor quality 

7 


Fragmentation 

0 


Total 

120 
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Fig. 2. Image sequence of an ice sphere containing 5% formic acid col¬ 
liding with a pure ice target (indicated) at 90° to the direction of travel at 
a velocity of 0.08 m s~ ] . The images were captured using beam splitter 
optics with the view on the left separated from the view shown on the 
right by 48.8°. The images are taken 4/107 s apart. 

Table 2. A break down of analysed collisions in terms of chemical com¬ 
position. The angles refer to the angle between the target and the direc¬ 
tion of travel of the particles. The ice spheres were composed of pure 
hexagonal crystalline water ice, 5% methanol in water and 5% formic 
acid in water. 


Collision type 

Pure water 

5% 

methanol 

5% 

formic acid 

Total 

Target at 90° 

7 

10 

18 

35 

Target at 60° 

8 

0 

0 

8 

No target 

2 

8 

5 

15 




Total 

58 


3. Analysis methodology 

The data analysis was performed as follows. The particles were 
manually tracked frame by frame in both views of the colli¬ 
sion; the co-ordinates from the two views were combined us¬ 
ing a transformation algorithm to give orthogonal (x, y, z) co¬ 
ordinates. Linear fits in each dimension yielded particle trajec¬ 
tories from which velocities before and after the collision and 
hence coefficients of restitution were calculated. The velocities 
were resolved into components normal and tangential to the col¬ 
liding surfaces to give normal and tangential coefficients of resti¬ 
tution. For the particle-particle collisions, normalised impact pa¬ 
rameters (b/R) were also calculated from particle trajectories. 
Where errors are shown, they were calculated using standard 
methods for the propagation of errors. 


4. Results and discussion 

4. 1. Coefficient of restitution 


Bouncing was the only collisional outcome observed in this 
study, corroborating previous results in Paper I and also those 
of HeiBelmann et al. ( 2010f . If ice particles have a critical ve¬ 
locity for the onset of bouncing, it will be below 0.01 m s 1 (the 
lowest velocity collision in the current study) for particles of this 
size. Fig.[3]shows the coefficients of restitution (s) as a function 
of relative impact velocity for this study compared the previous 
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Fig. 3. Coefficient of restitution as a function of impact velocity for 
collisions of 1.5 cm (diameter) ice spheres with a target at 90° to the 
direction of travel ( b/R = 0, shown by crosses), a target at 60° to the 
direction of travel ( b/R = 0.5, shown by squares) and no target (both this 
data (shown by circles) and data previously reported by HeiBelmann 
|et al.| l |20T0) (shown by stars)). The chemical composition of the spheres 
are indicated with colours (black for water, red for methanol, blue for 
formic acid and grey for the previous data). The rectangles (black for 
pure water, red for 5% methanol in water and blue for 5% formic acid 
in water) show the data encompassed by one standard deviation of the 
mean. 


results of jHeiBelmann et al.l ( |2010j > which were for collisions of 
pure water ice spheres of identical size. Apart from two outliers 
around 0.01 m s -1 , the data in the current study lies virtually 
within the range of the previous data. The two distinct veloc¬ 
ity groupings in our data are because of the different relative 
velocities obtained by colliding a particle with a stationary tar¬ 
get (the lower velocity region) and colliding a particle with an¬ 
other moving particle (the higher velocity region). The velocity 
range is similar for both studies (0.05-0.19 m s _1 for the cur¬ 
rent study (disregarding the two outliers around 0.01 m s -1 )) 
and 0.06-0.22 m s 1 for the previous study). The values of s 
are spread between 0.08 and 0.81 for our data and between 0.02 
and 0.84 for the previous data; the range is virtually identical 
for the current and previous results. There is no apparent differ¬ 
ence between the collisional properties of the different chemical 
compositions. The similarity of the target data with the previ¬ 
ous particle-particle collision data shows that colliding a sphere 
with a larger body has the same effect as colliding it with another 
sphere, a result which is not surprising considering the similar¬ 
ities in contact area during collision in both cases. Finally, the 
similarity of these particle-particle collisions with the previous 
work reiterates the fact that the brief opening of the chamber 
during the flight to remove the target did not affect the results; if 
significant frosting of the particle surfaces had occurred during 
this time, the coefficients of restitution would be lower (Hatzes 
let al.|19881 >. 

Where there was sufficient data, rectangles have been plot¬ 
ted to show how many points lay within one standard deviation 
of the mean (ignoring any outliers). These rectangles have been 
plotted independently for each experimental setup. For the target 
at 90° and for particle-particle collisions, there is a good overlap 
between the rectangles, demonstrating that there is no statistical 
difference between the datasets. We therefore conclude that the 
presence of methanol and formic acid in the ice at a level of 5% 
does not change the coefficient of restitution. For the target at 
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60°, technical difficulties prevented the capture of useful colli¬ 
sions for any other material than pure water, meaning that the 
effect of chemical composition cannot be tested for this flight. 

At this point it is instructive to consider why the inclusion 
of methanol and formic acid does not noticeably affect the col- 
lisional properties of the ice spheres. As mentioned in the in¬ 
troduction, methanol frosts are much stickier than water frosts 
between 110 and 150 K at atmospheric pressure (Bridges et al. 


1996]); this is most likely because of surface wetting near the 


melting point. However, surface wetting is not possible at the 
low pressures within our chamber for either methanol or formic 
acid. It is also useful to consider where the methanol is located 
within our ice spheres. Considering the methanol-water phase 
diagram ([Takaizumi & Wakabayashi (1997]), Deschamps et al. 


( 2010| >), at 255 K (the temperature of a standard kitchen freezer), 
the water will have frozen but the methanol will remain as a liq¬ 
uid. As the freezing process will freeze from the outside in, it is 
possible that all of the methanol is at the centre of the ice sphere 
due to exclusion from the water ice as it freezes, meaning that 
the surface is essentially pure water ice which would explain the 
similarity in collisional properties. Formic acid freezes at 282 
K, which is higher than the freezing point of water, although it 
is likely that the presence of the water will depress the freez¬ 
ing point of formic acid. At a concentration of 5%, we speculate 
that the formic acid will be distributed throughout the ice sphere 
due to the similarity in freezing points. Regardless of where in 
the spheres the methanol and formic acid is present, at a level of 
5%, the overall structure is likely to be dominated by crystalline 
water ice. It is likely for this reason that the collisional properties 
show no dependence on the composition of the ice. 

To determine the strength of any correlation between the 
variables, the linear Pearson correlation coefficient was com¬ 
puted: 


EEqCh - y)(xj - x) 

VE'LoO'/-j) 2 VE'Lo(-F-x ) 2 


(3) 


The number of points is given by n and the mean values of x 
and y are given by x and y respectively. This gives a measure of 
the strength of correlation between two variables, x and y, with a 
value of 1/-1 indicating perfect positive/negative correlation and 
a value of 0 indicating no correlation. The value returned is only 
significant where there is a significant number of data points and 
so only cases that meet this condition will be discussed. For the 
target at 90°, the data for pure water, methanol and formic acid 
shows no correlation between coefficient of restitution and im¬ 
pact velocity (values of 0.21, 0.06 and 0.05 respectively). The 
same is true for pure water with the target at 60° (correlation co¬ 
efficient of 0.13) and for methanol with no target (correlation co¬ 
efficient of 0.47); the rest of the datasets contain an insufficient 
number of points for any correlation to be statistically signifi¬ 
cant. These values indicate that the coefficient of restitution does 
not depend significantly on impact velocity which corroborates 


our previous work in Paper I and the work of HeiGelmann et al. 
( 20T0| . Here, as before, the surface roughness is likely to be the 
most important factor in explaining the range of coefficients of 
restitution. 

Fig.[4]shows the coefficients of restitution tangential and nor¬ 
mal to the colliding surfaces as a function of tangential and nor¬ 
mal impact velocity respectively. Calculating these gives further 
information about the distribution of energy after the collision. 
The normal coefficient of restitution gives information about the 
rebound of the particles and the tangential coefficient of restitu¬ 
tion gives information about particle scattering. The data for the 


target at 90° (b/R- 0) and the data for particle-particle collisions 
(b/R- 0.00-0.33) both have a small component of tangential im¬ 
pact velocity compared to the normal component (Fig. |4ji, b, e, 
f), which is to be expected for head on (or nearly head on) col¬ 
lisions. The tangential coefficients of restitution are spread from 
around 0.1 to 1.3 and from around 0.2 to 1.1 for particle-particle 
collisions. The spread for the normal coefficients of restitution 
is much less, from around 0 to 0.6 for the target at 90°, and 
from around 0 to 0.7 for particle-particle collisions. This shows 
that head on collisions can result in both scattering and rebound, 
with more energy going into scattering than rebound on average. 
For the target at 60° (b/R- 0.5), there is a much larger tangen¬ 
tial component to the impact velocity than for the target at 90°, 
although the normal component to the impact velocity is still 
larger (Fig. & and d). The spreads of tangential and normal co¬ 
efficients of restitution are very similar in this case, both around 
0.3-0.7. This demonstrates that for b/R- 0.5, the distribution of 
energy into rebound and scattering is roughly equal. In all cases, 
the chemical composition of the ice does not affect the distribu¬ 
tion of translational kinetic energy. 

Considering the highest value of s obtained in this study, a 
minimum of 29% of the particles’ translational kinetic energy 
(1 - s 2 with e = 0.84) is lost in the collision. From the image 
sequences of the collisions, it is clear that some of the energy 
is converted into rotational energy. A break down of the particle 
rotation is shown in Table [3] The overall values are given but 
dividing the spheres by chemical composition does not signif¬ 
icantly alter the proportions. Due to the spherical shape of the 
particles (and hence the lack of distinguishing marks), it was not 
possible to extract more quantitative rotational information. 

Table 3. The percentage of particles that rotate and do not rotate before 
and after the collision. 


Rotates (%) 

Does not rotate (%) 

Unclear (%) 

Before 4 

89 

7 

After 71 

10 

19 


The majority of the particles were fired from the pistons 
without any observable rotation (89%). The majority of the par¬ 
ticles did rotate after the collision (71%) but 10% did not. It is 
clear from this and previous results in Paper I that rotation can¬ 
not account for all of the particles’ energy loss. A computational 
study by |Zamankhan| ( |2010| > showed that most translational ki¬ 
netic energy in icy grain collisions is dissipated due to surface 
fracturing. We cannot confirm or refute this here because such 
surface fracturing would not be visible in our images. It is also 
possible that energy is converted into heat leading to desorption 
of surface material, which again is not detectable in our experi¬ 
ment and is beyond the scope of this paper. 

4.2. Effect of temperature 

Due to parabolic flight restrictions, it is not possible to cool the 
chamber with liquid nitrogen during flight. The flow of nitrogen 
must be stopped prior to take off, relying on the copper in the ex¬ 
periment to act as a heat sink. The chamber temperature slowly 
increases during the flight, giving us an opportunity to study the 
effect of temperature on coefficient of restitution. The data for 
the third flight of this campaign have been excluded due to the 
possibility of additional warming when the chamber was opened 
to remove the target. The temperature ranged from 131 to 160 K, 
increasing at a rate of 22.5 K li 1 for the first 20 minutes and at 
a rate of 7.0 K li for the remainder of the flight. Fig. [5] shows 
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Fig. 4. Coefficients of restitution normal and tangential to the colliding surfaces as a function of the normal and tangential impact velocity 
respectively. The chemical compositions are indicated by different colours and symbols: black plus symbols for pure water ice, red crosses for 
5% methanol in water, and blue circles for 5% formic acid in water, a. Tangential coefficient of restitution for the target at 90° to the direction of 
travel ( b/R = 0). b. Normal coefficient of restitution for the target at 90° to the direction of travel (b/R = 0). c. Tangential coefficient of restitution 
for the target at 60° to the direction of travel (b/R = 0.5). d. Normal coefficient of restitution for the target at 60° to the direction of travel (b/R = 
0.5). e. Tangential coefficient of restitution for particle-particle collisions (b/R = 0.00-0.33). f. Normal coefficient of restitution for particle-particle 
collisions (b/R = 0.00-0.33). 


the coefficient of restitution as a function of temperature. The 
data for the second flight (target at 60°) is at lower temperatures 
than the data for the first flight (target at 90°) because the avail¬ 
able collisions from this flight all took place earlier in the flight 
than those on the first flight. There is no correlation between the 
parameters (a correlation coefficient of -0.06) and hence we con¬ 
clude that temperature has no effect on coefficient of restitution 
over this temperature range. This corroborates our previous work 
in Paper I and again we conclude that this is due to the lack of 
surface melting at these temperatures and pressures. In order for 
surface melting to take place, higher temperatures and pressures 


would be required (around 250 K and a few mbar, considering 
the phase diagram of water ( Ehrenfreund et al7||2003) l). In this 
temperature and pressure regime, we would expect to see a re¬ 
duction in coefficient of restitution with increasing temperature 
owing to greater sticking forces at higher temperatures, as seen 
by Supulver et al. (T997]). 


5. Astrophysical implications 

It appears that inclusion of methanol and formic acid at a level 
of 5% has no impact on the collisional properties. However, 
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methanol freeze out in planetary rings and protoplanetary disks 
will happen at lower temperatures than those at which our freez¬ 
ing process took place, which could lead to a truly mixed ice. In 
this scenario, methanol may have a greater effect on collisional 
properties, particularly if it is present in higher concentrations 
near the ice surface. However, formic acid has not been detected 
above the 5% level and is likely to be evenly distributed through¬ 
out our ice sample, leading us to conclude that the inclusion of 
formic acid will not affect the collisional properties of the ices. 
Therefore, it seems unlikely that the inclusion of other species 
on the level of 5% will affect the collisional outcomes of icy par¬ 
ticles in ring systems or protoplanetary disks. 

Both the pure ice spheres and the ice spheres containing 5% 
methanol or formic acid yielded a broad range of coefficients 
of restitution. This corroborates previous results by HeiGelmann 
|et al.| ( |20T0| l and Paper I. Historically, the coefficients of restitu¬ 
tion determined by Bridges et al. (19841 have been most com¬ 
monly employed in numerical simulations of ring system dy¬ 
namics. However, our results show that since both elastic and in¬ 
elastic collisions take place concurrently, it is more appropriate 
to use a broad range of coefficients of restitution, evenly spread 
between 0.08 and 0.81. This conclusion is valid for both pure ice 
spheres and ice spheres containing 5% methanol or formic acid. 

Although surface melting as a sticking mechanism is not pos¬ 
sible at these temperatures and pressures for crystalline ice, a 
similar surface restructuring mechanism is likely to be possi¬ 
ble for amorphous solid water, which has been detected in star 
forming regions ( Smith & Wright|2011[ > and around young stel¬ 
lar objects (Scheg erer & Wolf||20 1 0| > , and is speculated to be 
present in protoplanetary disks. Amorphous solid water under¬ 
goes a glass transition to cubic ice at 137 K ( [Smith & Kay|1999] l 
which is within the temperature range of our experiments. Colli¬ 
sions at this temperature would be more likely to lead to sticking 
as bonds could be formed between the particle surfaces as they 
undergo restructuring and this has already been proposed as a 
sticking mechanism (|Supulver et al.|1997jl. 


6. Conclusions 

The main conclusions of our work are as follows: 


1. Sticking does not occur as a collisional outcome. Universal 
bouncing was observed, demonstrating that the critical ve¬ 
locity for the onset of bouncing of ice particles of this size 
(1.5 cm diameter) is less than 0.01 m s . 

2. No difference was observed between the coefficients of resti¬ 
tution for pure crystalline water ice, water with 5% methanol 
and water with 5% formic acid. The presence of methanol 
and formic acid in the ice does not affect the collisional 
behaviour in our experiments. We postulate that this is be¬ 
cause all samples have surface structures that are dominated 
by crystalline water ice. As in previous studies (Paper I and 
|HeiGelmann et al.| ( |2010) >), the surface roughness is the dom¬ 
inant feature of the particles, not their chemical composition. 

3. A broad range of coefficients of restitution was found with 
no correlation between this and impact velocity, replicating 
the results of previous studies (Paper I and 
( [2010] )). As before, this is thought to be 
anisotropic surfaces of the ice particles. 

4. At least 29% of the particles’ initial translational kinetic en¬ 
ergy is lost in the collision. Some of this energy is converted 
into rotational energy but this cannot account for all of the 
energy loss. 

5. Temperature did not affect coefficient of restitution over the 
range measured (131 to 160 K). 
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